DIAGONALIZATION OF AN INTEGRABLE DISCRETIZATION 
OF THE REPULSIVE DELTA BOSE GAS ON THE CIRCLE 
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Abstract. We introduce an intcgrable lattice discretization of the quantum 
system of n bosonic particles on a ring interacting pairwise via repulsive delta 
potentials. The corresponding (finite-dimensional) spectral problem of the 
integrable lattice model is solved by means of the Bethe Ansatz method. The 
resulting eigenfunctions turn out to be given by specializations of the Hall- 
Littlewood polynomials. In the continuum limit the solution of the repulsive 
delta Bose gas due to Lieb and Liniger is recovered, including the orthogonality 
of the Bethe wave functions first proved by Dorlas (extending previous work 
of C.N. Yang and C.P. Yang). 
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1. Introduction 



The non-ideal Bose gas with deha-potential interactions is a system of n one- 
dimensional bosonic particles characterized by a Hamiltonian given by the formal 
Schrodinger operator 

H^-^ + g Y, S{xj~Xk). (1.1) 

l<j^k<n 

Here xi , . . . , a:„ represent the position variables, A := d'^_^ -I- • • • -I- d^^ , d refers to the 
delta distribution, and g denotes a coupling parameter determining the strength of 
the interaction. For g > the interaction between the particles is repulsive and 
for g < it is attractive, whereas for 5 = the model degenerates to an ideal 
one-dimensional boson gas without interaction between the particles. 

The eigenvalue problem for the above Schrodinger operator with periodic bound- 
ary conditions — i.e. for particles moving along a circle — was solved by Lieb and 
Liniger by means of the Bethe Ansatz method |LLj . The corresponding spectral 
problem for particles moving along the whole real line was considered subsequently 
by McGuire |Mcj . Since the appearance of these two pioneering papers, the exactly 
solvable quantum models under consideration have been the subject of numerous 
studies; for an overview of the vast literature and an extensive bibliography we 
refer the reader to Refs. jMal IG4I IKBII lAKI |S]. Further generahzations of the 
one-dimensional quantum n-particle system with delta-potential interactions can 
be found in Refs. |G3IIGullHO[IDi| . where analogous quantum eigenvalue problems 
are studied in which the permutation-symmetry is traded for an invariance with 
respect to the action of more general reflection groups [Bl IHu| . and also in Refs. 
[AFKLIXkIIGGIIHLPJ . where the delta-potential interaction between the particles is 
replaced by more general zero-range point-like interactions (involving combinations 
of S and S' type potentials) |A-H| . 

An important (and notoriously hard) problem connected with the Bethe Ansatz 
method is the question of demonstrating the completeness of the Bethe wave func- 
tions in a Hilbert space context. For the non-ideal Bose gas on the line with a 
repulsive delta-potential interaction {g > 0), the spectrum of the Schrodinger op- 
erator is purely continuous. The completeness of the Bethe wave functions was 
proved for this case by Gaudin |G1I IG2| . For the corresponding system in the at- 
tractive regime (g < 0), the completeness problem is much harder as multi-particle 
binding may occur thus giving rise to mixed continuous-discrete spectrum. In this 
more complex situation the completeness of the Bethe wave functions was shown 
by Oxford [O], with the aid of techniques developed by Babbitt and Thomas in 
their treatment of an analogous spectral problem for the one-dimensional infinite 
isotropic Heisenberg spin chain [Tl IBT| . When passing from particles on the line 
to particles on the circle the nature of the system changes drastically, as the con- 
finement to a compact region forces the spectrum of the Schrodinger operator to 
become purely discrete. In this situation the completeness of the Bethe wave func- 
tions was proved for the repulsive regime by Dorlas jDoj . In Dorlas' approach the 
question of the completeness is first reduced to that of the orthogonality of the 
Bethe wave functions. This orthogonality is then shown to hold with the aid of 
quantum inverse scattering theory jKBIj . combined with previous results of C.N. 
Yang and C.P. Yang pertaining to the solution of the associated algebraic system 
of Bethe equations (determining the spectrum of the Schrodinger operator under 
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consideration) lYYj . For the attractive regime such progress has yet to be made: 
the question of the construction of a complete eigenbasis for the non-ideal Bose gas 
on the circle with delta-potential interaction remains (to date) open. 

By exploiting the translational invariance and the permutation symmetry the 
eigenvalue problem characterized by the Hamiltonian H (|l.l|l reduces — in the case 
of bosonic particles moving along a circle of unit circumference — to that of the free 
Laplacian 

-AV- = E^j, (1.2) 

acting on a domain of wave functions 'ip{xi, . . . , a;„) with support inside the alcove 

A = {x e M" I xi H h a;„ = 0, xi> X2> ■ ■ ■ > Xn, 2:1 - x„ < 1}, (1.3) 

and subject to normal linear homogeneous boundary conditions at the walls of the 
alcove of the form 

(a,^. -a,^.+i-5)V'L^_,^.^^^o = 0, J = 1,..., 71-1, (1.4a) 

(a.„-a,, -5)^l,^_,^^i = o. (1.4b) 

(The parameter E represents the energy eigenvalue.) The purpose of the present 
paper is to study a discretization of the eigenvalue problem in Eqs. (|1.2|l - l|1.4b|l . 
Throughout we will restrict attention to the repulsive parameter regime 5 > 0. 

More specifically, we study the eigenvalue problem for an integrable system of 
discrete Laplacians acting on functions supported on a regular lattice over the al- 
cove A (|1.3|l . and subject to repulsive reflection relations at the boundary of the 
lattice. Since the alcove is compact, the lattice in question is finite; hence, our 
discrete eigenvalue problem is finite-dimensional. We solve the eigenvalue problem 
at issue by means of the Bethe Ansatz method. The resulting Bethe eigenfunc- 
tions turn out to be given by specializations of the Hall-Littlewood polynomials 
|M2[ IM3| . The orthogonality and completeness of these Bethe eigenf unctions arises 
as an immediate consequence of the integrability (which permits removing possible 
degeneracies in the spectrum of the Laplacians). As a byproduct, the Lieb-Liniger 
type Bethe eigenfunctions for the eigenvalue problem in Eqs. H1.2|l - ljl.4b|l are re- 
covered via a continuum limit. The orthogonality of the latter eigenfunctions (and 
thus eventually — because of Dorlas' results |Doj — also the completeness) are in our 
approach immediately inherited from the corresponding orthogonality results for 
our discretized lattice model. In this connection it is probably helpful to recall that 
the original proof of the orthogonality due to Dorlas |Doj also involves a discretiza- 
tion, which arises however in a fundamentally different way from the one employed 
here. In a nutshell: Dorlas arrives at the orthogonality through a continuum limit 
of the (second) quantization of the Lattice Nonlinear Schrodinger Equation intro- 
duced by Izergin and Korepin |KBI| . whereas here — in contrast — we study a rather 
more elementary quantum lattice model characterized by a direct discretization of 
the Schrodinger operator in Eqs. (|1.2(l - l|1.4b|l itself. 

The paper is structured as follows. In Section |21 the discretization of the eigen- 
value problem in Eqs. H1.2|l - ljl.4b|l is formulated. In SectionOthe eigenfunctions are 
constructed by means of the Bethe Ansatz method. The associated Bethe equations 
are solved in Section^and the orthogonality and completeness of the corresponding 
Bethe wave functions is demonstrated in Sectional Finally, the continuum limit is 
analyzed in Sectional 
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2. Discrete Laplacians on the Alcove 

In this section we introduce a system of discrete Laplacians on a finite lattice 
over (a dilated version of) the alcove A (|1.3|l . For this purpose it will be convenient 
to borrow concepts and notation from the theory of root systems. Here we will 
only need to deal with the simplest type of root systems: those of type A. For 
further background material concerning root systems the reader is referred to Refs. 
iRllH^. 



2.1. Preliminaries. Let ei,...,e„ denote the standard basis of unit vectors in 
R" and let (•, •) be the (usual) inner product with respect to which the standard 
basis is orthonormal. The alcove A (|1.3(l constitutes a convex polyhedron in the 
center-of-mass plane 

^ = {xeM" I (x,e) =0}, e = ei + --- + e„, (2.1) 

which is bounded by the n hypcrplanes 

(2.2a) 

= l,...,n-l, (2.2b) 



Eo = {xeE 


1 (x,q;o) = 1}, 


Ej ^{xeE 


|(x,a,)=0}, 


where 




ctQ := ei — e„ and 


«j - e, - e,_ 


Specifically, we have that 




A = {xeE\ {x, ao) 


<1, (^■a,)2 



j^l,...,n-l. (2.3) 

j = l,...,n-l}. (2.4) 

The vertices (i.e. corners) of the polyhedron A are determined by the intersections 
of all choices of n — 1(= dim{E)) out of the n hyperplanes Eg, . . . ,En-i- These 
vertices are given explicitly by the origin and the vectors 

a;j:=eiH hej-^(eiH he„), j = l,...,n-l. (2.5) 

Indeed, the vectors Wi, . . . , a;„_i all lie on the hyperplane Eq (|2.2a|) and constitute 
a basis of £^ that is dual to the basis ai, . . . , ocn-i (in the sense that (wj, cxk) = Sj^k, 
where Sj^k denotes the Kronecker delta symbol). 

Let ro : E ^ E he the orthogonal reflection in the hyperplane Eq (|2.2a() and let 
rj : E ^ E, j = 1, . . . , n — 1 be the orthogonal reflections in the hyperplanes Ej 
(I2.2b|) . The action of these reflections on an arbitrary vector x E E is oi the form 

ro(x) = x+ (1 - (x, Q:o))ao, (2.6a) 

rj(x) = X — (x, q:j)q:j, j = l,...,n— 1. (2.6b) 

From these two formulas it is readily inferred that for j e {l,...,n— 1} the re- 
flection rj swaps the j^^ and (j + 1)*'' coordinates of x, and that tq swaps the 
flrst and the last coordinates followed by a translation over the vector ao- Hence, 
the reflections ri, . . . ,r„_i generate an action of the permutation group 5„ on E, 
and the reflections rp, . . . , ?'n-i generate and action of the afflne permutation group 
Sn = Sn <>< Q, which is the semidirect product of the permutation group 5„ and the 
lattice of translations 

Q:=Spa.ii^{ai,...,otn-i). (2.7) 

For any (afflne) permutation a E Sn, one deflnes its length £(a) as the minimal 
number of reflections needed for decomposing a (non-uniquely) in terms of the 
generators: 

o- = ^ji'"i2---^i. (2-8) 
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(where ji, . . . ,ji e {0, 1, . . . , n — 1} and with the convention that the length of the 
identity element is equal to zero). 

The polyhedron A (|2.4|l constitutes a fundamental domain for the action of 5„ 
on E. More specifically, for each x G E the orbit iS„(x) intersects A precisely once. 
Let us denote by CTx G 5n the unique shortest affine permutation such that 

ax(x) e A. (2.9) 

Let us fix a positive integer m. Below it will often be convenient to employ a 
dilated version of the polyhedron A*-™-* := mA rather than A itself. Throughout 
we shall distinguish by means of a superscript (m) the corresponding boundary 
planes, boundary reflections, and the elements of the affine permutation group 
Sn := Sn X (mQ) C 5„ generated by these reflections. That is to say, the 
dilated alcove A^™' is bounded by the hyperplanes E"^' := rnEj, j = 0, . . . , n— 1; 
the orthogonal reflections in these hyperplanes act as rj" (x) :— ■mrj{x/rn), j = 
0, ...,n— 1. (So E^™' = Ej and r^^ — rj if j > 0.) The affine permutation 
(Tx G Sn mapping a vector x. E E into the fundamental domain A^"^' is given 
by the action cr^ {^) '■— 'tiCTx /,„. (x/m) . 

2.2. Laplacians. The orthogonal projection of Z" C M" onto the center-of-mass 
plane E ()2.1|l is given by the lattice dual to Q (|2.7|l : 

V := {XeE\yaeQ: {\,a) eZ} (2.10a) 

= Spanz(a;i,...,u;„_i). (2.10b) 

It is clear from Eq. H2.1()a(l that Q is contained as a sublattice in V. Hence, the 
action of the affine permutation group 5„ in E maps the lattice V into itself (cf. Eqs. 
I)2.6a|l . (j2.6b|) '). The intersection of the lattice V with the dilated polyhedron A*^™^ 
provides a finite grid T-''-™-' over A^"^' containing its vertices and muJi, . . . , mujn-i'- 

P(") := {Ae-P I Ae A*")} (2.11) 

= {kiuji H kn^iujn^i I fci, . . . , K-i e Z>o, fci H h fc„_i < m}. 

We are now in the position to define a system of n — 1 Laplace operators acting in 
the space C{T'^™'^) of complex functions V-" '■ T^^™-* -^ C 

Definition (Laplace Operator). To each basis vector u}^. from Eq. (|2.5|) we asso- 
ciate a corresponding Laplace operator D^' : C(P*^'"-') -^ C{V^"^'') defined by its 
action on an arbitrary function i/j : T't™' — > C of the form 

(4"V)a:- Y. V^^+- (2.12a) 

with the boundary convention that ior fi £ V\ T''™) 

V'M-^''"'^^^"'V<.«(^), (2.12b) 

where t denotes a real (coupling) parameter (and the length function f ("*)(•) refers 
to the minimal number of reflections in the decomposition of an affine permutation 
in Sn in terms of r^ , ■ ■ ■ , ?'„_i)- 
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Roughly speaking, the value of L^ ip i^i ^ point A e 7^(™) is equal to the sum of 
the values of tp in all neighboring points of the form A + i', where v runs through 
the orbit of w^ with respect to the action of the permutation group 5„. When \ + v 
lies outside p'™' the value of V'A+f is governed by the boundary convention in Eq. 
(I2.12bl) . It is instructive to clarify the nature of this boundary convention somewhat 
more in detail by decomposing Ox+i, in terms of the elementary reflections in the 
hyperplanes bounding A^"^' . 

Proposition 2.1 (Boundary Reflection Relations). Let A € T't™). The boundary 
convention in Eq. I|2.12b|l amounts to the requirement that \tv € Sn{<^k) for which 

\ + v e ■p\7'('") 

tV'.l)(,^.) ^f{X + u,c,)<0 0->0), ('-'^"^ 

or equivalently 

, = f i^A+i^-ao */ {\ "o) = m and {v, ao) = 1 

[#A+i/+a, if (A, Oij) = and {u, cxj) = -1 (j > 0). 

Proof. Let A e -pf") and u e 5„(a;fc) such that \ + u <E V \ V^'^\ Then there 
exist i £ {0, . . . ,n — 1} such that (A + i/, a.j) is either > to if j = or < if 

j > 0. Geometrically, this means that the hyperplane E "^ separates X + v from 

^x+vi^ + I/) G 7'(™). Let us write fi :— rj™ {X + u). Then for any such j we have 
that 

We will now use induction on the length of crx+ij ^'~' prove that the boundary 
convention in Eq. (|2.12b|) and the boundary reflection relation in Eq. (|2.13a|l are 
equivalent. Indeed, upon assuming l)2.12bp it is clear that 









which amounts to (|2.13a() . Reversely, upon assuming H2.13a|l and invoking the 
induction hypothesis we see that 

which amounts to Ij2.12b|l . To finish the proof of the proposition it remains to check 
that the boundary reflection relations in Eqs. H2.13a(l and l|2.13b|l are equivalent. 
For this purpose it suffices to notice that the requirements that A £ J^i™-) and 
u £ Sn{^k) imply that < (A, ckj) < m and that —1 < {v^olj) < 1 for j — 
0, . . . , 71 — 1. Hence (A + v, ao) > m iff (A, cxq) = m and {v, a.^) = 1, in which case 
Tq -^ (A + I/) = A + 1/ — CKo, and furthermore for j > one has that (A + u, ctj) < 
iff (A, ctj) — and (v, ctj) = —1, in which case rJ-™^(A + i^) = X + i^ + ctj . D 

It is clear from the proposition that the boundary convention in Eq. (|2.12b|) 
amounts to a normal linear boundary condition at the hyperplanes bounding A^"^' . 
The coupling parameter t determines the nature of this boundary condition; for 
|i| > 1 the boundary term (i.e. the interaction between the particles) is attractive 
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whereas for |i| < 1 it is repulsive. For t = Q and t ~ 1 we are dealing with Dirichlet 
type and Neumann type boundary conditions, respectively. 

Applying the boundary convention to those contributions in the sum over the 
translated orbit A + 5„(u;fe) corresponding to lattice points outside the grid 7'(™) 
gives rise to a closed formula for the action of the Laplace operator in which the 
value of L^ ip at the point A is expressed completely in terms of the values of 
■0 at the neighboring points of the form X -\- u £ 7'(™). To make this procedure 
completely explicit we shall need some further notation. Let R be the orbit of the 
basis cci, . . . , q:„_i with respect to the action of the permutation group 5„ and let 
R"*" be the part of the orbit that expands nonnegatively with respect to this basis: 

R = {gj - Gfe I 1 < jV fc < n}, R+ = {gj - Gfe I 1 < j < fc < n}. (2.14) 

Proposition 2.2 (Explicit Action of the Laplace Operator). The action of the 
Laplace operator Lj^ (k £ {1, . . . , n—1}) on an arbitrary grid function i/j : V^"^' -^ 
C is of the form 

iLt^^h= E ^i"^'/'A+., (2.15a) 






whe 



yW:= TT ^ \ ^ n ^-— ^^, (2.15b) 

c>;GR+ aeR+ 

(A,a)=0 (A,c>;)=r?j. 

(i/,ct) — 1 (t/jCt) — — 1 

with p := I]„eR+ "^/^ = ^^i H V t^n-i- 

Proof. It follows from (the proof of) Proposition 12.11 that for A G V^™^ and v G 
Sn{i^k) such that X + u e 7'(™) the coefficient V^"^' of ipx+u in (L^" V)a is of the 

form "^a- i^ ^'^' I where the sum is over those affine permutations a in Sn that are 
generated by iterated application of reflections of the type figuring in Eqs. (|2.13a|l . 
(|2.13b|) . These affine permutations are precisely the a € Sn for which a{X) — X 
and cr(A + v) ^ pi™-) ^ or equivalently, a[X + i>) ^ X + u. In other words, the 
coefficient is equal to the Poincare series of the quotient of the stabilizer subgroup 
^nx '■— {'^ ^ '^'" I '^W — ^} ^-nd the stabilizer subgroup 5^"^' n S^^x+v^ i-^- 



Viz'- y. ^""^-^ - ^ ^^""^ ...... . (2.16) 



^o-GSi^'n^i'"'. 



It follows from a general formula for the Poincare series of (affine) Weyl groups due 
to Macdonald IMl| (cf. Corollaries (2.5) and (3.4)) that the Poincare series in the 
numerator and denominator of Eq. (|2.16() admit product representations given by 

y /™'M^ TT ^-^ / / TT 1^^ 1_ (2.17) 

o-GS*"\' aGR+ aeR+ 

"^'■^ {X,oc)=0 {X,a)=m 
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and 

y /-H'^)^ TT ^"^ TT Lll (2.18) 

"■^ ".A+L. {A,a)=0 (A,a)=m 

(A+i^,q;)=0 (A+j/,a)=m 

respectively, which — upon inserting in Eq. (|2.16(l — gives rise to Eq. (|2.15b|) . D 



Remark. The stabiHzer subgroups iS„'x ^'^'^ ^nX ^ ^nX+u ^"^ ^^^ proof of Proposi- 
tion |^| consist of (direct products of) permutation groups. It is well-known (and 
readily seen by induction) that the Poincare series of the permutation group S^ 
admits the product representation J|-j^(l — t-') /(I — t) (cf. e.g. Ref. [M2], Chapter 
III, §1). With the aid of this latter product formula it is not so difficult to verify 
Eqs. (|2.17|l . (|2.18|l (and thus Eq. (|2.15b(l 'l directly (i.e. without invoking the much 
more general results of Macdonald |M1| cited in the proof) . (We thank the referee 
for making this point.) 

2.3. Hilbert Space Structure. We shall now endow the function space C(7'''™^) 
with a inner product, turning it into a (finite-dimensional) Hilbert space Ti*-™' := 
£2(-p(m)^^(m)) characterized by a positive weight function A^™) : "P^") -> (0,oo). 
To this end it will always be assumed from here onwards that the coupling parameter 
t lies in the repulsive regime 

-l<t<l (2.19) 

(unless explicitly stated otherwise). For two arbitrary functions ip^cj) & C(P'-'"') the 
inner product in question is then defined as 

(V',0)('"^= J2 V^A^Af\ (2.20a) 

AGP'"' 

where the weight function is given by 

, . T 1 - /(P-°) -,-r 1 — f^-iP'") 

4"^- n hk^ n hh^i^- ('-'^^^ 

aeR+ aeR+ 

(A,a)=0 {X,a)=m 

(Notice in this connection that the restriction of the coupling parameter to the 
repulsive regime (|2.19|l ensures that the values of the weight function A^^ are 
indeed positive for all A G T^C™).) 

Proposition 2.3 (Adjoint). The Laplace operators L]J^ and L^_j. (k £ {1, . . . , n— 
1}) are each others adjoints in 7i(™) 

VV',0eC(7'(™)) : (LfV^^)^™^ = (V,4-U)^'"^- (2-21) 
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Proof. The proof hinges on the exphcit formula for the action of LjT^ in Proposition 



12.21 Elementary manipulations reveal that 



A 

AGPt™) 

= E E vs[:^VA+.^A(r) 

A+i^e-p<") 

= E E ^.^i"l^;:^A(r) 



where we have used: (i) the substitution \ ~ fi ^ f, (ii) the identity 

^i"LAirl = C-lAM, (2.22) 

and ('iizj the facts that the coefficient Kj"_j^ is real and —ujk G Sni(^n-k)- To infer 
the identity in Eq. H2.22|l one observes that for /x £ T-"*-"*^ and v € ^^(wfc) such 
that fi — uG •p('") both sides V^l,,^ ^ll-i, and V^^^, A)I"' reduce — upon cancehng 
common terms from the numerator and denominator — to 

li 1 - ll + {p,a) li J^ _ ^l+n-(p,a) ' 
Q;eR+ aeR+ 

(/x,a)— {fJ,,a)—7n 

O 

3. Bethe Ansatz Eigenfunctions 

In this section the eigenfunctions of the Laplacian L™ (|2.12a() . Ij2.12b|l are 
constructed by means of the Bethe Ansatz method of Lieb and Liniger jLO IMal 

in4llEET| . 

3.1. Bethe Ansatz. If we ignore boundary effects for a moment and interpret 
-^fc H'^-f^a|l (without the boundary convention (|2.12b|) ') as a Laplacian acting on 
functions ip : V ^ C, then clearly the plane wave ipxi^) = exp(z(A, ^)) with wave 
number ^ G E/{2t:Q) constitutes an eigenfunction corresponding to the eigenvalue 
Ek{$,) = J2ves (u> ) 6xp(*('^!^))- The Bethe Ansatz method aims to construct 

the eigenfunctions tp : 7^^™' -^ C for the operator L^ (|2.12a|) with the boundary 
convention l|2.12b|l via a suitable linear combinations of plane waves (corresponding 
to the same eigenvalue Ek{£,))- This prompts us to look for eigenfunctions given 
by a linear combination of plane waves exp(j(A, $,„)), a G 5„ with coefficients such 
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that the boundary conditions in Proposition 12.11 are satisfied. Specifically, we will 
employ a 5„-invariant (in ^) Bethe Ansatz wave function of the form 

*>^(^) = 7?7T E i-irC{Cy<f+^^^^\ i e 2rfnt(A), (3.1a) 

where |, := a(0, (-1)" := clct(a) = (-l)^('^), 

and Int(A) = {^ G £? | (|, ccq) < 1, (|, ccj) > 0, j = 1, . . . , n - 1}. (The condition 
that ^ e 27rInt(A) guarantees that the denominator d{^) is nonzero.) 

Proposition 3.1 (Bethe Wave Function). The Bethe Ansatz wave function "^xi^) 
()3.1a|l . (|3.1b|) satisfies the boundary reflection relations in Provosition \2.1\ for i = 
1 , . . . , ?i — 1 provided that 

C(0= n (l-te-'<°'«>) (3.2) 

aeR+ 

(or a scalar multiple thereof). 

Proof. Let A e 7^(™) and u e Sn{(^k) such that (A + t^,OLj) = —1 for some j G 
{1, . . . , n — 1}. Equating 

*A+.(0 = ^ E (-l)^C(C)e'^^+^+'^'«'> 

to 

= ^ E (-l)^C(C)e*<°-«">e'^^''+^+'^'«"> 

leads to the relation 

Y, {-~irC{i.r)^t E (-irC(C.)e'<"-«--> Vre5„. 

Because rj stabilizes p + \ + u (i.e. rj G 5„_p4-A+i/), the latter relation can be 
rewritten as 

^ (-l)-[C(C.)-C(r,(C,))] 

""^Sn.p + A + L- 

= t J2 (-l)'^[C(l..)e*<°-«"> -C(r,(0)e-'<«-«^^>]. 

a-i(a.)eR.+ 
By varying A and u it is seen that this relation implies that 

C{^)-C{r,m ^ t[C(Oe*<°-«> -C(r,(0)e-<°-«>] 

(as an identity in ^) for all reflections r-,- , j = 1, . . . , n — 1, or equivalenty (assuming 
C(^) is nontrivial in the sense that it does not vanish identically) 

C(0 _ l-te-'<«^'«> 



C{rj{^)) l-ie*(°^-.«> 



j = l,...,n- 1. 
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Hence C(^) must be of the form 

aeR+ 

where Co(^) denotes an arbitrary 5,i-invariant overall factor (i.e. cq{^^) = co(^), 

Vcr e Sn). □ 

3.2. Bethe Equations. By pulling the overall factor 6{^) inside the sum and 
exploiting the anti-invariance 6{^^) — {—1)°'S{^), the Bethe wave function ^a(C) 
(I3.1a|l . Ij3.1b|l . with coefficients C(^) taken from Eq. H3.2II . passes over to 

From this expression it is clear that — for A fixed — the Bethe wave function amounts 
to a Hall-Littlewood polynomial in the spectral parameter ^ M2, M3 . We will 
now derive conditions on the spectral parameter such that the Bethe wave function 
satisfies the boundary reflection relations in Eqs. 12.13a|l . (|2.13b|) for j — 0. 

Proposition 3.2 (Bethe System). The Bethe Ansatz wave function ^a(C) I|3.3(I 
satisfies the boundary reflection relations in Proriosition \2.1\ for i = if the spectral 
parameter ^ G 2TTlnt{A) solves the algebraic system 

^^raip^i) = (L^ ) TT L^ , v/3eR. (3.4) 

ctGR 

{a,/3) = l 

Proof Let us write C(^) := Y[^^^+ \1*JL7(1°o and let A G T''™) and v e 5„(a;fc) 
such that {\ + v,ao) = rn + I. Equating 

creSn 

to 

yields the relation 

E (l-te~''<°°'«->)C(|,)e'<^+"'«-> =0, 
cres„ 
or equivalently 






(1 - te-'^<°-'-o'"(«')>)C(^r(C))e'^^^+'^"""'<«-") = 0, 



where Tq e Sn denotes the orthogonal reflection Tq (x) = x — (x, Q:o)ao- The 
latter equation translates to 



a-\ao)eR+ 



J2 ((l-te-<«-«">)C(0 + 

«o)eR+ 

(1 _ ie^(«o,«„))c(r^o)^^^^^g_(™+iH„o,«j^g^(A+.,cj ^ 0^ 
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which is satisfied if 

(1 - ie-*<«°'«">)C(^„) + (1 - te*<°°'«->)C(r^°^(|„))e-("+i)*<"°'«-> = 
for all (T G 5„, or equivalently (assuming C{^^) =/= 0) 

^(™+i).(/3,,) _ _C(rf(0) 2^|e^^ 

- C(C) l-ie-^'^^^^' "' 

where /3 := ^^^(ao). The proposition now follows upon inserting 

{a,ao)>Q {oc,l3)>a 



4. Solution of the Bethe Equations 



D 



In this section the Bethe System in Proposition 18 . 21 is solved using a variational 
technique due to C.N. Yang and C.P. Yang JYYllM^inii rkBT . 

4.1. Solution. The following theorem provides (the existence of) a sequence ^„ G 
27rlnt(^) of solutions to the Bethe system in Proposition 13.21 labeled by vectors 
(playing the role of quantum numbers) /i. G 7^(™). 

Theorem 4.1 (Bethe Vectors). For each jj, G T''™) there exists a (unique) Bethe 
vector ^^ G 2TrInt{A) = {x G .E | (x, cto) < 27r, (x, a.j) > 0, j = 1, . . . , n — 1} such 
that ^^ satisfies the system in Eq. (|3.4|) . Moreover, these Bethe vectors have the 
following properties: 

(i) $/x' = iti «/ O'f^d only if fj.' = fi, 

(ii) ^ depends smoothly on the boundary parameter t G (—1, 1), 
(''') ^t.^7^iP + fJ')fort = 0. 

4.2. Proof. In standard coordinates the Bethe system of Proposition 13.21 reads 



'-(?.-?.)= TT ^ ~ '^'^ 1 TT '-'^ 1 (4 1) 



i<e<n i<e<n 



for 1 < j ^ k < n. This overdetermined system of n(n — 1) equations in the 
variables ^i , . . . , ^n is equivalent to the system of n equations 

^'"'^-^^ n e-(0-^.)-, ' J-^^---^-, (4.2) 

l<i<n 

where c ^ denotes an overall constant factor that we can scale to 1 by means of 
the translation £^j —^ ^j — im~^ log c, j = 1, . . . n. Picking thus c = 1 and taking the 
logarithm of both sides recasts Eq. 14. 2|) in the additive form 

Kj+ J2 0i^j~^e)^2Trmj, j = l,...,n, (4.3) 

l<l<n 

i^3 
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where m = (mi, . . . , m„) £ Z" and 

0{x) := {1-t^) I {l-2tcos{x)+t^)-^dx (4.4a) 

Jo 

f \ +t X \ 

= 2arctan tan(-) | (4.4b) 

Here the branches of the arctangent function and those of the logarithmic function 
are to be chosen in such a way that (i) 9{x) 14.4bll . 14.4c() is quasi-periodic: 9{x + 
27r) = 9{x) + 2tt, and (ii) 9{x) (|4.4b|l . (|4.4c|l varies from —n to tt as a; varies from 
— TT to TT (which corresponds to the principal branch). We notice that this choice of 
the branches ensures that 9{x) (|4.4b|) . (|4.4c|l is smooth on the whole real axis and 
strictly monotonously increasing. 

Lemma 4.2. For each n-tuple m = (?7ii, . . . , to„) G Z", there exists a unique 
vector ^(m) = (^i(m), . . . , ^„(m)) solving the system in Eq. H4.3I) (with 6{x) of the 
form in Eqs. H4.4a|l - H4.4cll ). Furthermore, this solution ^(m) depends smoothly on 
the boundary parameter f G (— 1, 1). 

Proof. Let 

n -. n n 

nei,...,U:=f E^'+2 EQ(«^-^^-)-2^E"^^-^^' (4.5) 

where 9(a:) :— J^ 9{x)Ax. Clearly the solution(s) of the system in Eq. (|4.3I) coincide 
with the critical point(s) of the (smooth) function F(^i, . . . , <^„). The Hesse matrix 
of V is given by 

H,.u = ^1^ = ("^ + E ^'(0 - 6) I '5,,fe - 9'{^, - a-), 1 < J, fc < n, 




where 9'{x) = (1 — t^)(l — 2icos(a;) + i^) ^ > 0. It is readily seen that this Hesse 
matrix is positive definite: 

n n ^ n n 

E Hj^kXjXk =™E^J + 2 ^ ^'i^J~^k)ixj -Xkf >m'^x]>Q 

j,k=l j=l j,k=l j = l 

(for any nonzero vector x e M"). The function y(^i, . . . , ^„) is thus strictly convex, 
i.e., it admits at most one critical point: a global minimum. That such global 
minimum ^(m) indeed exists in our case is immediate from the observation that 
V{S,i, ■ ■ ■ ,S,n) — > +0O when ||^|| — > oo. (Notice in this connection that Q{x) — > +00 
for X -^ ±00.) We thus conclude that the system in Eq. (|4.3I) has a unique solution 
^(m) (given by the global minimum of V). It remains to check that the position of 
this global minimum depends smoothly on the boundary parameter t. To this end 
we notice that the integrand of 9{x) (|4.4a() (which, incidentally, coincides with the 
generating function for the Chebyshev polynomials) is analytic in i for |i| < 1, and 
thus so are the function V{£,i, . . . ,£,n) and the system of equations for the global 
minimum in Eq. H4.3|l . The smoothness in ^1 , . . . , ^„ and t — combined with the fact 
that the Hessian det{Hj,k) is positive (i.e. nonvanishing) — now guarantees that the 
solution ^(m) to the latter system must be smooth in t € (—1, 1) by the implicit 
function theorem. D 
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The following lemma shows that the ordering between the components of the 
solution ^(m) coincides with the ordering of the components of the labeling vector 
m. 

Lemma 4.3. Let m e Z" and let ^(m) be the associated solution of Eq. (|4.3|) 
detailed in Lem,m,a \4.S\ Then for rrij > ra^ the following inequalities hold 

2Tr(niA — nil.) , , , , 2'K(mA — mi.) 

^-^ -^ < m - ^ m < ^-^ -^ , 4.6a 

where K±{t) := ^j^^^isa > 0. So, one has in particular that 

mj > TOfe =4> ^j(m) > ^fc(m) and m^ = mj. =4> ^j (m) = ^^ (m) . (4.6b) 

Proof. Let m in Z" with ruj > nik. Subtracting the fc*'* equation from the j*'* 
equation of the system in Eq. 14.3|l yields that 

n 

™fe- - 6) + E(^(^J- - ^^) - ^(^fc - ^^)) = 27r(m, - m,). (4.7) 

Since the r.h.s. of this identity is nonnegative and 9{x) is strictly monotonously 
increasing, it follows that ^j(m) > Cfc(m)- Furthermore, from the formula d{x) — 
e{y) = (1 - t'^) /^(l - 2icos(a;) + t'^)-^dx (cf. Eq. (|Oa|l 'l it is immediate that 

K,-\.{t){x — y) < 9{x) ~ 9{y) < n^{t){x ~ y) for x>y. 

Application of this upper and lower bound so as to estimate the terms in the sums 
on the l.h.s. of Eq. 14.7|l . now gives rise to the inequalities 

(m + nK+{t)){^j{m) - Cfe(m)) < 2'n{mj ~ ruk) < (m + nK_(i))(^j(m) - 6(m)), 

which completes the proof of Eq. 14.6a(l (and thus also that of Eq. Ij4.6b|l ). D 

The next lemma improves the upper bound on the distance between the Ci(m) 
and ■^fc(m) stemming from Lemma |4. 31 in the situation that the distance between 
rrij and mk is smaller than n + m. 

Lemma 4.4. Let m G Z" such that mj^rrik < n+m and let ^(m) be the associated 
solution of Eq. H4.3(l detailed in Lemma \4.S\ Then 

ej(m)-a(m)<27r. (4.8) 

Proof. Subtracting the k*^ equation from the j*'' equation of the system in Eq. 
(14.31) leads — upon recalhng that 9{x) is odd — to (cf. Eq. H4.7|) ') 

n 

™(e, - Ck) + ^(^(ej - Ci) + 0{it - 6)) - 27r(m, - mfc). (4.9) 

1=1 
If ^j — ^k > 27r, then the average of ^j — ^g and ^i — f fc is > tt. Hence 9{^j — S^g) + 
9{^g — ^fc) > 27r, in view of the fact that 9(x) is strictly monotonously increasing 
and 0(7r + a;) + 0(7r — a;) = 27r. Plugging this estimate in Eq. (|4.9|l reveals that 
^j(m) — Cfc(m) > Stt implies that 2TT{mj — mk) > TO(^_,(m) — ^^(m)) + 27rn > 
27: {m + n), which completes the proof. D 

We will now piece the results of Lemmas 14.2114.41 together, so as to arrive at a 
proof for Theorem 14. II 
Let m G Z" such that 

mi > m2 > ■ ■ ■ > rUn and mi — m„ < m + n, (4.10a) 



E g. l4^ 
E g. l4'3l 
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and let ^(m) be the associated solution of Eq. (|4.3(l detailed in Lemma 14.21 It 
follows from Lemmas 14.31 and 14.41 that 

a(m)>e2(m) >--->e„(m) and ^(m) - 6i(m) < 2^. (4.10b) 

Let us define 

/x:=m (m,e)e-p, ^^ := |(m) (|(m),e)e, (4.11) 

where (recall) e = ei + • • • + e„ and p = Ui + ■ ■ ■ + uin-i- In other words, /x is the 
orthogonal projection of m onto the ccnter-of-mass hyperplane E (12.11) translated 
by — p and ^ is the orthogonal projection of ^(m) onto E. The inequalities in 
Eqs. (|4.1()a(l and l|4.10b|l ensure that /x e T''™) (|??TT|l and that |^ G 27rlnt(74) 
(cf. Eq. (|2.4|l ). respectively. It is furthermore clear that by varying m we can 
reach any lattice point /x e T''™). Indeed, for /x — kiuii + • • • + fc„_ia;„^i with 
kj e Z>o and fci + • ■ • + fc„_i < to we may pick the components of m equal 
to nij — (/I + p, aj + ■ • ■ + OLn-x) — kj + ■ ■ ■ + fc„_i + n — j, j = 1, . . . , n. It 
is not difficult to check that the assignment p ^ ^ is indeed well-defined (i.e. 
/i, = p,' => ^ = ^ ,) and one-to-one (i.e. ^ = ^ , => p, = fi'). Indeed, one has 
that 

|(m) - l(m') e 27rZe 

m m' e Ze 
^=> p = fi'. 

Since it is obvious that ^ inherits from ^(m) the smooth dependence on the 
boundary parameter t and the property that its components solve the system in 
Eq. (|4.2|l (and thus the Bethe system in Eq. H4.1|l l. this proves Theorem 14. II up to 
Property (ii). It remains to check Property (Hi), which states that for t = the 
Bethe vectors are given by £„ = „^„ (p + m), A* G pi"^) , To this end we simply 
observe that Lemma [4.31 implies that for t = 

27r 

0(m) - Cfe(m) = — ■ (toj - TOfc), 

n + m 

whence the statement follows by varying m subject to the constraints in Eq. H4.10a|) 
and projecting onto the center-of-mass plane with the aid of Eq. H4.11(l . 

5. DiAGONALIZATION 

In this section we will combine the results of Sect ions |2}S to arrive at an orthog- 
onal basis for the Hilbert space 7i(™) = £^(7'(™), A^™)), consisting of a complete 
set of joint eigcnfunctions for the (commuting) Laplace operators i^™ , . . . , i„™i- 

5.1. Spectrum and Eigenfunctions. The following theorem provides the eigcn- 
functions of our Laplace operators in terms of Hall-Littlcwood polynomials special- 
ized at the Bethe vectors ^^, p S p('"). 

Theorem 5.1 (Spectrum and Eigenfunctions). For special values of the spectral 
parameter, given by the Bethe vectors ^^, fi G 'P*-'"' in Theorem \4-l[ the Bethe 
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wave function ^a(^) I|»^.'^|I constitutes an eigenfunction of the Laplace operator Lj!^^ 



(|2.12a|l . (|2.12b|) . i.e. for any k € {1, . . . ,n~ 1} 
where the eigenvalue is of the form 



4"^ *(!;.) = i?fc(C^)*(U' (5-la) 



Ekii)^ E exp(z(iv,0) (5.1b) 

rand ^{^^) ^ 0). 

Proof. Clearly the Hall-Littlewood polynomials ^a(€) I|3.3|) satisfy the identity 
J2ues„{ujt)^><+'^(i) ^ -Efc(|)*x(l) (because ah of the plane waves tpxi^a) = 
exp(i(A,^^)), a e 5„ do so). Moreover, since the specialized Hall-Littlewood poly- 
nomials ^E'A(^^l), fJ' G T'f™) also satisfy the boundary convention in Eq. (j2.12b|) 
in view of Propositions 12.11 13.11 13.21 and Theorem 14.11 the stated eigenvalue equa- 
tion follows. It remains to check that ^a(^„) does not vanish identically. For this 
purpose it is enough to observe that for A = 0: 



i — te ""' (i) ^ «(^) (m) -n- i-t 



*o(o = E n \^u:u - E ^''^ - n 



1 - t<p'«) 



where we have used (i) a rational function identity and (ii) a product formula for 
the Poincare series of the permutation group that arc both due to Macdonald jMl] 
(cf. Theorem (2.8) and Corollary (2.5), respectively). It is clear from the product 
formula on the r.h.s. that ^o(€) > for — 1 < i < 1, whence ^a(€^) indeed 
constitutes a true (i.e. nonzero) eigenfunction in Ti'™). D 

5.2. Orthogonality and Completeness. Theorem 15 . II provides as many eigen- 
functions as the dimension of the Hilbert space (indeed, dim(7i(™)) = :^'p("^) — 
\ m~ ))■ "^^"^ following theorem confirms our expectation that these eigenfunc- 
tions actually form an orthogonal basis for the Hilbert space in question. Alterna- 
tively, one may think of this theorem as describing a novel system of discrete (dual) 
orthogonality relations for the Hall-Littlewood polynomials. 

Theorem 5.2 (Orthogonality and Completeness). The Bethe wave functions 

*(l;x), MGT'^™) (5.2a) 

constitute an orthogonal basis of 'H^"'-> : 

V/x,/.'e7'(™): (*(^^),*(V))^'"' = |° V^^'' (5-2^) 

>0 i//x = /i'. 



Proof. Since Lj^ and i„_j, are each others adjoints in 7i(™) by Proposition 12.31 

it is clear that (i^^^(l^), *(|^,)>^™^ = (^(I^J^^I^^^m'))^"^' ^^ applying 
Theorem 15.11 and using the fact that Ek{^) — En-k{^), this equality is readily 
rewritten in the form 

{E,{^^) - i?fc(V))(*(|^),vI^(|^0>'™' = 0. (5.3) 

Theorem 14 . 1 1 now guarantees that for fi y^ /j,' the associated Bethe vectors ^ and 
^„/ are distinct in 27rInt(A). Moreover, since the elementary symmetric polynomials 
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Ei{^), . . . , En-i{^) separate the points of 27rInt(A) (as they generate the fuU alge- 
bra of trigonometric polynomials on 2ttA spanned by the 5„-invariant Fourier basis 
S«eS„(A) 6xp(i(/i,|)), A e V), this imphes that in this situation i?/j (^^ ) ^ £'fc(^^,) 
for a certain value oi k E {1, . . . ,n— I}. We thus conclude from Eq. 1)5. 3|) that the 
inner product (^(^„), *(^„'))*-™'' must vanish if /i 7^ /i'. Finally, for fi ^ fi' the 
inner product yields the squared norm of the Bethe wave function ^(^„) in Ti''"^\ 
which is positive as ^(^^) 7^ by (the proof of) Theorem 15. II D 

5.3. Integrability. From the previous results it is seen that our Laplace operators 
model a finite-dimensional quantum system that is integrable in the following sense. 

Theorem 5.3 (Integrabihty). The Laplacians l["'\ ..., L,^™\ (|2.12a|) . (j2.12bp con- 
stitute n — 1 (~ dim{E)) mutually commuting operators in the Hilbert space Ti.^"^'. 
Furthermore, any operator L : Ti,^™-' -^ Ti'™-* that commutes with all of the Lapla- 
cians Lj^ , . . . , i„™i lies in the polynomial algebra C[Lj^™ , . . . , L„™j^]. 

Proof. The commutativity of L" , . . . , L^li is immediate from the fact that the 
operators are simultaneously diagonalized by the basis ^(^„), M G 7^^™^ of Ti^™) (cf. 
Theorems l5 . 1 1 and l5 . 2)1 . The property that any operator L : TiS™-^ -^ Ti^™' that com- 
mutes with Li \ . . . , L„_]^ is necessarily algebraically dependent of L^ \ . . . , L„_i 
hinges on the fact that the eigenvalues i?i(|^), . . . , £'„-i(^^) separate the elements 
of the eigenbasis ^(£„), M S p('") (cf. also the proof of Theorem 15.2(1 . Indeed, it 
is immediate from this that L is diagonahzed by ^(^^), ^J■ G T^^™). In other words, 
that there exist a function E^ '■ {€/i}(j,e75(m) ^ C such that 

L^{i^.) = EL{i^)-^{i^). VMe7'(™). (5.4a) 

Since the Bethe functions ^(^^), A* G T'^™) form an orthogonal basis of Ti^'"^ 
we have (by transposition) that the Hall-Littlewood polynomials ^a(^), A e 7'(™) 
form a basis for the space of complex functions on the spectral set {£,u} u,^-p{m) upon 
specialization. In particular, there exist (unique) complex coefficients cx, A S 7-"^™^ 
such that 

EL{i^)= E ca*a(|^), V/xg7'("). (5.4b) 

Aep('") 
Furthermore, from the well-known property that the elementary symmetric poly- 
nomials ii^i(^), . . . , £'„_i(^) (|5.1bp generate the space of symmetric polynomials it 
is clear that there exist a polynomial P^ S C[i?i, . . . , i?„-i] such that 



^ ca*a(0-^l(^i(0,---,^«-i(0)- (5.4c) 



E 

ASP'""' 

It follows from Eqs. (|5.4a() - H5.4c|l and Theorem 15.11 that the operators L and 
Pl{L\^' , . . . ,L]^\) coincide on the basis ^(^^), A* G p^."^) , Hence, we conclude 
that L = PUl["'\ . . . , Lt\) e C[l["'\ ..., Lt\]. □ 

The Laplace operators L]^ , . . . , L^i are not self-adjoint in general in view 
of Proposition 12.31 As a consequence, the spectrum in Theorem 15.11 is generally 
complex- valued. Within the commuting algebra C[Ll^' , . . . , L„_i] there exist how- 
ever many operators that are self-adjoint. For example, the alternative generators 
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-(™) -'- /rC"') _L r(™) 



L'n,l ■■- ^{Lr+Lr-i)^ fce{l,...,[n/2]}, (5.5a) 

4? - ^(4"^-4"l), ke{l,...,[{n-l)/2]}, (5.5b) 
are self-adjoint and have real spectrum of the form 

ErA^^) = E cos((i.,|^)), /xeT'^™), (5.6a) 

EiAi^) = E sin((«.,^^)), /xep(™), (5.6b) 



respectively. The real subalgebra R[L]^j', . . . , ij^r„/2ii -^/"i j • ■ • i -^/7(n-i)2l] '^°^^" 
sists of all operators L : Ti^'") — s- Ti^™) such that (^ij L commutes with all of the 
Laplacians L^"'' , . . . , L„™\ and (ii) L is self-adjoint. 

One of the simplest positive operators in this real subalgebra is given by 

ij('") :==nId-L^"'. (5.7) 

In standard coordinates the explicit action of this operator on an arbitrary wave 
function f/; G Ti^™' is of the form (cf. Proposition 12. 2(1 

(i/(™)V)A = nV-A - i E ^TA^-A+.^-i ^ V-^^x.^^, (5.8a) 



where 



^^A 



^.:a 



2 

l<j/<n l<j/<n 

X+UjeV'-'"' A-i^j6-p<"' 



n -T^ n ^T^FT^' (5-8b) 



j<k<n l<k<j 

\k=\j Afc=A,+m 



11 1 _ p-k 11 l-tn+J-k ' (5-^^) 



l<fc<j j<k<n 

Afc — A. Ai. — A,— 7^ 



and i/j = ej — (ei -I- • • • -|-e„)/ri, j = 1, . . . , n (so fi , . . . , !/„ consist of the orthogonal 
projection of the standard basis ei, . . . . e„ onto the center-of-mass plane E (|2.1|l ). 
The spectrum of _ff(™) is built of positive eigenvalues E{^^), fi e 7^(™) with 

n 

E{^) = J2{l^cos{^j)). (5.8d) 

i=i 

The operator ij(™) (|5.8al) - H5.8c|l serves as the Hamiltonian of our lattice n-particle 
model. Below we will verify that in a continuum limit this lattice Hamiltonian 
tends formally to the Hamiltonian of the n-particle delta Bose gas on the circle. 

6. Continuum Limit 

In this final section we first review the solution of the eigenvalue problem for 
the Laplacian in Eq. (|1.2|l . with wave functions supported inside the alcove A 
(|1.3I) subject to repulsive boundary conditions of the form in Eqs. (|1.4al) . (|1.4bp 
(i.e. with g > 0). Our formulation amounts to the center-of-mass reduction of the 
seminal results due to Lieb and Liniger |LL| (Bethe wave functions), C.N. Yang and 



DISCRETE DELTA BOSE GAS ON THE CIRCLE 19 



C.P. Yang ..YY (Bcthc vectors), and Dorlas |Doj (orthogonality and completeness). 
Next we will show how this solution of the eigenvalue problem for the Laplacian in 
the alcove can be recovered from the corresponding solution of our discrete lattice 
model via a continuum limit. 



6.1. Eigenfunctions. In the notation of Section |21 the eigenvalue problem in Eqs. 
ifT^ - lfTIbl reads 

-AV' = Eip, -xG A, (6.1a) 

with 

((VxV',ao)+5V')|xeEo=0, (6.1b) 

((VxV',a,>-#)|xe£, -0, j = l,...,n-l (6.1c) 

(where Vx refers to the gradient). Let us define 

C := {^eE\{ta,)>0, j^l,...,n-l}, (6.2a) 

V^"°^ := {kiUJi + --- + kn-iuJn-i\ku...,k„^ieZ>o}. (6.2b) 

Theorem 6.1 (Bethe Wave Functions lEEj)- The Bethe wave function* 
with the spectral parameter $, £ C (|6.2a|) solving the Bethe system 

{a,/3> = l 

constitutes a solution to the eigenvalue problem in Eqs. (|6.1a(l " (l6.1c|l corresponding 
to the eigenvalue E = E'^'^\() := (|,|). 

It is instructive to recall briefly the essence of the proof of Lieb and Liniger in 
the present notation. Firstly, it is clear that the linear combination of plane waves 
^(°°^(x, ^) constitutes an eigenfunction of —A with eigenvalue E'^°°\$^). It remains 
to check that the boundary conditions are also satisfied. The boundary condition 



*This explicit form of the expressions for the coefficients of the Bethe wave function is due to 
Gaudin l;G.3 G4 . 
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in Eq. H6.1c|l is inferred by the following computation for x e Ej: 



- » E ( n ^il47^)'"-=-' 

ae5„ aeR+ ^ ' ^°"' 

where in Steps (i) and (ii) one exploits that ncteR+ lal )^ ^^"^ {^ji^a) ^''"^ 

symmetric and skew-symmetric, respectively, with respect to the action of Tj on 
^^, combined with the symmetry (x,rj(^g.)) — (x,^^.) (since rj(x) = x if x G £^j). 
Finally, the boundary condition in Eq. I|6.1b|) requires that for x G £^o 

Manipulations similar to those in the proof of Proposition 13.21 reveal that this 
relation holds when the spectral parameter solves the Bethe system in Eq. Ij6.3b|l . 

Theorem 6.2 (Bethe Vectors |^). Let g > 0. For each /x G p(°°) (|OH)l there 
exists a (unique) Bethe vector £„ G C (|6.2a(l such that ^ satisfies the system in 
Eq. (|6.3b|) . Moreover, these Bethe vectors have the following properties: 

(ii) ^ depends smoothly on the boundary parameter g > 0, 
(Hi) £ ^ 2n{p + /x) for g -^ +oo. 



In standard coordinates the Bethe system in Eq. I)6.3b|l reads 

^ -11 zg-6+C^ -^^ ig-6 + a' ^ ^ ^ 

l<£<n -^ ^-^ ^ l<£<n '^ ^ ^ 

for 1 < J 7^ fc < n, or equivalently (upon exploiting the translational invariance) 



,*«. = rr ^g + ^j - c^ 



---U'Bt^' , = l,...,n. (6.5) 
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In the additive form the latter system becomes 

0+ E 0^^\^,-i,) = 2iTm,, j = l,...,n, (6.6) 

l<l<n 

with m = (jTii, • ■ ■ , mn) € Z" and 

0(°°)(x) = 2g f\x^+g^y'dx (6.7a) 







= 2arctan(-) (6.7b) 

9 

= .logf^. (6.7c) 

\ig-xj 

It was shown in |YYj that for any m S Z" the Bethe system in Eqs. H6.6ll - H6.7c|l 
has a unique solution ^(m) given by the unique global minimum of the strictly 
convex function 

^ n ^ n n 

y(-)(Ci,...,en):-2E^I+2 E e^°^H^,-ik)-2nJ2m,^,, (6.8) 

j=l j,k=l 3 = 1 

with Q^°°\x) := J^ 6^°°'{x)dx. By projecting the solutions ^(m), corresponding to 
vectors m G Z" with mi > TO2 > • • • > m„, orthogonally onto the center-of-mass 
plane E the statements of Theorem 16. 21 readilv follow (cf. also Section^. 



Theorem 6.3 (Orthogonality and Completeness jDop . The Bethe wave functions 
^(°°''(x, ^ ), /x G 'p(°°) form an orthogonal basis for the Hilbert space T-C'^°°^ :— 
L'^(A, dx.) (with inner product {(f),il;)^°°' := J.(j){x.)ip{x.)dx.), i.e. 

V/x,/x'G7'(-): {^^°'H^^)M^\i^')y^'' = l\ Y^^': (6.9a) 

y>Q if n^ ^l' 

and 

(0, *^°°^ {^^)) = 0, V/x e p(°°) =^0 = 0. (6.9b) 

Below we will infer that this center-of-mass reduction of Dorlas' orthogonality 
relations can be recovered via a continuum limit from the corresponding results 
pertaining to the discrete lattice model in Section |S1 It was moreover shown by 
Dorlas that the orthogonality of the Bethe wave functions for the repulsive delta 
Boson gas implies their completeness jDol Section 3] . In other words, the complete- 
ness in Theorem 16 . 31 follows from the orthogonality (upon a cosmetic adaptation of 
Dorlas' arguments to our center-of-mass situation). 

6.2. Orthogonality. In order to perform the continuum limit let us from now on 
rescale the coupling parameter t putting 

t = e-9/", g>0. (6.10) 

For any x in (the closure of) C (|6.2a|l we define an integral approximation [x] e 
V''°°^ (lOKll of the form 



[x] := [(x, ai)](jJi H h [(x, a„_i)]u;„_i, (6.11) 

where [x] denotes the integral part of x e M>o obtained through truncation. With 
these notations we are in the position to embed the Hilbert space of lattice functions 
^M = f('p'^"^)^ A^"^'>) into L'^{C, dx)hy means oi&lineaviniectionJ^"'^ : 7^(™) ^ 
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L^(C,dx) that associates to a lattice function tf> : 'p(™) ^ C a staircase function 
J(") (0) : C -> C of the form 

(J(™)0)(x) := I ^^^[-1 ^^' ["^^] ^ ^*"^' (6.12) 



for [mx] ^P(™). 

It is not difficuh to see that the staircase function j('") (0) has support on a bounded 
domain inside the dilated alcove (1 + —)A. This support shrinks towards (a subset 
of) A for TO — > oo. It is also not difhcult to deduce from this definition that 

/'(j('"V)(x)(J('")V')(x)dx = c„,™ J2 -^aV^^™'. (6-13) 

where c„,™ = Vol(a;i, . . . ,a;„_i)/TO"~i = l/{m"'~^y/^). Let ^(™)(x,|) be the 
staircase embedding of the Hall-Littlewood polynomial ^a(^) (|3.3|) 

^('"nx,!) := (J("^1'(0)(x) (6.14) 



A™^,*; ,(0, 



-9/™ p-''(a,€„> 



The following lemma states that, for to — > 00, the rescaled staircase function 
^j,(m)^jj. j_^-j l|(j ]^4II converges pointwise to the Lieb-Liniger Bethe wave function 

^(°°'(x, ^) H6.3a|l when x lies in the interior of A and to zero when x lies outside 
A. 

Lemma 6.4. For any ^ e C, one has that 

™— 00 10 ifx€C\A. 

Proof. The lemma readily follows from the explicit expression of the staircase 
wave function on the third line of Eq. 16.14|l . together with the observation that 

lim,„^oo ^[ZL] = 1 if X e Int(A) and lim,„^oo ^[Zl.] = if x G C\ A, and the fact 
that hm^^oo ;^["^x] = x. D 

Let us fix a /x G 'pi°°) and pick m sufficiently large so as ensure that fi £ p(™). 
We denote by ^[j" ^^'^ ^iT" the associated Bethe vectors detailed in Theorem 14. II 
and Theorem 16. 21 respectively. 

Lemma 6.5. For any fi e 1^^°°', one has that 

hm TO|(:")=^(r^. (6.16) 

Proof. Let mj := {fi, a.j) + • • • + {fi, ctn-i) + n — j, j = 1, . . . ,n. Then ^J^"' and 
^iT" correspond to (the projections onto the center-of-mass plane of) the (unique) 
global minima of V{^i, ...,^n) ^M and V'-°^\^i, ...,^n) JHIHI), respectively. The 
rescaled Bethe vector to.^^^ thus corresponds to the global minimum of the func- 
tion l^'™''(^i, . . . ,^„) := TO.y(^i/TO, . . . ,^„/to). The lemma now follows from the 
observation that for to —> 00 the strictly convex function ^^'"^(^i, . . . ,^„) tends 
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to V^°^'{£,i, . . . ,^n) uniformly on compacts (which imphes in particular that the 
global minimum of the y('™) converges to the global minimum of V'^°°^). D 

The proof of the orthogonality in Theorem 16.31 now hinges on the following 
proposition. 

Proposition 6.6. For all fi, fi' £ 'P'^°^\ one has that 

lim / *(")(x,|("^)*^'"Hx,€L"^)*f (6.17) 

m^ca Jc 

vi/(°°)(x,|(r')*(^H^a[?^^. 

Proof. It is clear from (the proof of) Lemma 16.41 and from Lemma 16.51 that the 
integrand and support of the integral on the l.h.s. converges pointwise to the inte- 
grand and support of the integral on the r.h.s. To see that the integrals themselves 
converge accordingly we write 

j ^(™)(x, !)*(") (x,|')dx 

where C(|) = n„eR+ '-\l^.Ct'' ■ After substituting | := ^^T^ and |' := ^J'^ 
the proposition follows for ttt, — > oo upon invoking Lemma 16.51 and the dominated 
convergence theorem of Lebesgue. Indeed, one has that 

and 

^(m) ^ fl if xe Int(A) 

["^1 ' [0 ifxeC\A 

pointwise for m ^ oo, and that |e'^<[""l^""(«M'"')-"-'(C'')>| = 1, |a["^,| < 1. D 



Proposition 16.61 can be rephrased as 



lim / (j('")vI/(|M))(x)(j(™)*(^^™'))(x)dx. 

The r.h.s. of this limiting relation vanishes when fi ^ fi' in view of Eq. H6.13|l and 
Theorem 15.21 whence the orthogonality in Theorem 16 . 31 follows . 

6.3. Hamiltonian. We will now wrap up by verifying briefly that formally the 
Hamiltonian i/(™) (I5.8a|) . (|5.8c|) converges in the continuum limit to the Hamil- 
tonian of the repulsive delta Bose gas on the circle. It is quite plausible that with 
a somewhat more in-depth analysis in the spirit of Ref. one would be able to 
show that this convergence of the Hamiltonian is in fact in the strong resolvent 
sense, but we will not attempt to do so here. 

Let H^°°' be the self-adjoint extension in 7^^°°^ = ^^(A, dx) of the Laplace 
operator —A with boundary conditions of the from in Eqs. (j6.1b|) . (|6.1cf) . and let 
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H^™-* be the following rescaled staircase embedding of the operator 7?'^™^ H5.8a|l - 
(lO^ ini2(c,dx): 

h(") = 2m2j(™^if(™)(j(™))^in('"\ (6.18) 

where II^™) : L^(C,dx) -^ L^(C, dx) denotes the orthogonal projection onto the 
finite-dimensional subspace of staircase functions J'^™^(7i'^™^) C L^(C,dx). It is 
clear that 

jj(oc)^(oo)(^(oo)^ ^ ^'°°Hl^°°^)*'°°^(l|r^), (6.19a) 

with £;(°°)(|) = (1,1), and that 

hM*M(^M) =^('")(4'"^)*^'"n4'"^), (6.19b) 

where £'(™)(|) := 2m'^E{^) with ^(|) given by Eq. ifOdjI . From Lemmas lOI 
andEElit follows that lim^^oo *'™Hx,^^™^) = *(°°Hx,l[r^) pointwise for x e 
Int(A) and that lim,„^oo ^'"'H^IT'^) = ^'""H^IT^)- In other words, for m ^ oo 
the eigenfunctions, the eigenvalues, and the eigenvalue equation for H^™^ in Eq. 
I|6.19b(l converge pointwise to the eigenfunctions, the eigenvalues, and the eigenvalue 
equation for h'""-* in Eq. (|6.19a(l . respectively. 

Acknowledgments. Thanks are due to M. Bustamante and to S.N.M. Ruijse- 
naars for several helpful discussions. 
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